Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 February 2008 



(MN MfeK style file vl.4) 



Refining the Oort and Galactic constants 



Rob P. Oiling 1 *, Michael R. Merrifield^ 

1 Dept. of Physics and Astronomy, 

University of Southampton, Southampton S 017 IB J, U.K. 



Accepted 1998 February 2 



ABSTRACT 



The local stellar kinematics of the Milky Way offer a useful tool for studying the 
rotation curve of the Galaxy. These kinematics - usually parameterized by the Oort 
constants A and B - depend on the local gradient of the rotation curve as well as its 
absolute value (0o), and the Sun's distance to the Galactic center (Ro). The density of 
interstellar gas in the Milky Way is shown to vary non-monotonically with radius, and 
so contributes significantly to the local gradient of the rotation curve. We have there- 
fore calculated mass models for the Milky Way that include this component, and have 
derived the corresponding radial variation in the Oort constants. Between 0.9i?o and 
1.2i?o the Oort functions A(R) and B(R) differ significantly from the general ~ Qq/R 
dependence. Various previously-inexplicable observations are shown to be consistent 
with these new predictions. For example, these models may explain the ~40% differ- 
ence between the values for 2ARq derived from radial velocity data originating in the 
inner and outer Galaxy ( Merrificld 1992). They also go some way toward explaining 
the different shapes of the velocity ellipsoids of giant and dwarf stars in the solar neigh- 
bourhood. However, a consistent picture only emerges if one adopts small values for 
the radius of the solar circle (Ro = 7.1 ± 0.4 kpc) and local circular speed (Oo = 184 
± 8 km s~ ). With these Galactic constants the Milky Way's rotation curve declines 
slowly in the outer Galaxy; V rot (20 kpc) = 166 km s _1 . Our low value for the distance 
to the Galactic center agrees well with the only direct determination of Ro (7.2 ± 0.7 
kpc, Reid 1993). Using these Galactic constants, we fi nd that the proper motion of Sgr 
A* is consistent with the observational constraints ( Backer fc Sramek 1987| ; Backer 



1996; Reid 1998). Simple analytic arguments as well as detailed calculations show that 



the radial velocities and proper moti ons of our best fit model are entirely consistent 
with the radial velocities of Cepheids ( Pont, Mayor fc Burki 1994 ) and the Hipparcos 
measurements of their proper motions ( Feast fc Whitelock 1997] ). 



Key words: Galaxy: structure - Galaxy: kinematics and dynamics - Galaxy: solar 
neighbourhood - Galaxy: fundamental parameters - Galaxy: stellar content - Cepheids 
- ISM: general 



1 INTRODUCTION 

Since the first measurements of the Milky Way's rotation 
almost seventy years ago, a great de al of effort has been ex- 
pended in studying its kine matics (e.g., Fich & Tremaine 



1991, and references therein). The simplest kinematic quan- 



tity that we can determine is the rotation curve, O(-R), which 

* E-mail: olling@astro.soton.ac.uk 
t E-mail: mm@astro.soton.ac.uk 



measures the speed of circular orbits as a function of radius. 
Unfortunately, our location within the Milky Way means 
that even this simple dynamical quantity is difficult to estab- 
lish. In particular, the modest uncertainties in our knowledge 
of the distance to the Galactic center (Ro = 7.7 ± 0.7 kpc; 
Reid 1993) and the rotation speed at the solar circle (Oo = 
200 ± 20 km s ; Sackett 1997) lead to significant uncertain- 
ties in the inferred form for 0(i?) (Fich & Tremaine 1991). 
Measurements of stellar kinematics provide one approach to 
constraining the form of the rotation curve. The quantities 
measured are parameterized by the Oort functions, 
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For a flat rotation curve we see that both A and B are 
inversely proportional to R, with smaller A and larger B 
values in the outer Galaxy. The quantity B is determined 
by studying the dependence of stellar proper motions on 



galactic] longitude, while A can be calculated trom analysis 
of either proper motions (fii) or radial velocities (v r ): 



neighbourhood. Unfortunately, the available observations of 
the Milky Way's rotational motion are not good enough to 
allow one to calculate the derivatives in Q(R) directly from 
the data. Instead, some suitable functional form must be fit- 
ted to the data, and the derivatives of the model can then 
be calculated for comparison with the Oort constants. 

The customary technique for fitting a model to rotation 
curve data for an external galaxy is to decompose the galaxy 
into different mass components: a stellar bulge and disk; an 
azimuthally-symmet ric molecular and atomic gas disk: and 
a dark matter halo ( ;.g,, van Albada fe Sancisi 19861 Ik enq 



A d sin 21 
A cos 21 + B 
4~74 



(3) 
(4) 



with d the distance to the object (in units of kpc), I Galactic 
longitude, A and B in km s -1 kpc" 1 , and hi in milli arcsec 



(mas) per year (Mihalas & Binney 1981). Note that the 



above equations for v r and fii are only correct for << 1. 

It is also possible to determine A — B = Q/R essentially 
independently from the individual values of A and B from 
proper motion surveys in the directions of I = 90 degrees 
and I = 270 degrees. In these directions, objects have a 
rather limited range in Galactocentric radius, and are thus 
only slightly affected by the radial dependence of the Oort 
functions. 

Finally, the combination — B/(A — B) can be estimated 
from the shape of the velocity ellipsoid of random stellar 
motions. Although the Oort functions vary with Galactic 
radius, stellar kinematic observations have only determined 
the Oort A and B constants, which are the values of these 
functions in the solar neighbourhood. Thus, in practice they 
tell us about the local shape of the rotation curve. However, 
it is important to note that the range over which the rotation 
curve parameters are effectively averaged depends on the 
range in Galactic radii over which the stellar observations 
are made. 

The orbits of stars in the solar neighbourhood can be de- 
scribed as an epicyclic motion due to their random motions 
with respect t he local stan dard of rest, superimposed on a 
circular orbit (Oort 1965). The length of the semi-major 
axis of the epicycle in the radial direction (a ep i) depends on 
the velocity in the radial direction (II): 



= 0.511/ -B(A - B), 



(•») 



which equals ~350, ~700, and ~1200 parsec for Cepheids, 
early-type stars, and late-type stars, respectively. Thus, the 
random motions of stars introduces some form of radial 
smearing of kinematical properties which precludes signif- 
icant variation of kinematical properties on scales smaller 
than a ep i- 

Perhaps the greatest difficulty in interpreting the esti- 
mated values for the Oort constants comes from the presence 
of the derivative terms in equations (|lj) and (^): the values 
for the Oort constants depend on any local bumps and wig- 
gles in the rotation curve as well as its global form. Any 
attempt to interpret the Oort constants must therefore con- 
sider the local topography of the rotation curve in the solar 



1987 ; Bcgcman 



1989 



Oiling 199,' 



Lake fe Feinswog 1989; Broeils 1992 



piling 1996a ). The free parameters in this 
mass decomposition are then constrained by requiring that 
the rotation curve for the model provides the best match 
possible to the observed kinematic data. In this paper, we 
seek to apply the same basic technique to the kinematic data 
for the Milky Way in order to produce a dynamical picture 
consistent with the observed Oort constants. 

Although the dominant contributors to the total mass 
come from the stellar components and the dark halo, the 
contribution from the gas cannot be neglected in this ap- 
plication. While the stars and dark matter are believed to 
be fairly smoothly distributed with radius, the distribution 
of molecular and atomic hydrogen (Eb and H I respectively) 
often show density enhancements like rings and arms. Such 
irregular distributions will produce a contribution to the to- 
tal rotation speed that can vary non-monotonically with ra- 
dius. Thus, although these components do not contribute 
much to the amplitude of the rotation curve, they can have 
a large effect on its local derivative, and hence the predicted 
values of the Oort constants. 

We start by presenting the Milky Way's 'observed' ro- 
tation curve and its uncertainties in §g| In section ^, we cal- 
culate the radial distribution of gas in the Milky Way, and 
show that it does, indeed, significantly affect the analysis 
of the Oort constants. Section ^ presents the mass compo- 
nent decomposition of the Milky Way's rotation curve, and 
shows that the resulting models are completely consistent 
with the constraints from the Oort constants as long as we 
adopt values for Ro and Oo at the lower end of the accept- 
able range. In section ^[ we discuss the implications of these 
results. And we summarize in section ^. 



2 THE MILKY WAY'S ROTATION CURVE 

We have constructed the rotation curv e of the inner G alaxy 
from the tangent point measurements (Malhotra 1995). For 



the outer Galaxy, we have used deter minations of the W(R ) 
curve [Eqn. (0)] based on H2 regions (|Brand fe Blitz 1993| ) 
and the thickness of the gas layer QMerrifield 1992| ). As 
mentioned in the introduction, the rotation curve inferred 
from radial velocity observations depends significantly on 
the adopted values of Oo and Ro, and so we constructed 
about 50 "observed" rotation curves for 5 choices of Ro (6.1, 
6.8, 7.1, 7.8, and 8.5 kpc) and about 10 values of Oo between 
170 and 230 kms -1 . All of these rotation curves are con- 
sistent with the observed radial velocity measurements in 
the inner and outer Galaxy. Generally speaking, small val- 
ues for Oo and/or _Ro yield declining or flat rotation curves, 
while the rotation curves constructed for large Oo and/or 
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Figure 1. Nine possible rotation curve (dots) for the Milky Way which are consistent with the tangent point data (Malhotra 1995), the 
radial velocities of H II regions (Brand & Blitz 1993) and the H I W(R) curve (Merrifield 1992). For clarity, we did not plot the error 
bars on the observed data points (a few km s , and ~ 15 km s _1 in the inner and outer Galaxy, respectively). Each panel corresponds 
to a different choice for Ro and So, where all Ro-©o combinations are consistent with the Kerr & Lynden-Bell values at the 1.3-cr level. 
From left to right the Galactocentric distance Ro increases from 7.1 through 7.8 to 8.5 kpc. In the bottom row we present our best-fitting 
model (left panel; R = 7.1 kpc, 6 = 184 km s -1 , A-B = 25.92 km s _1 kpc _1 ; see §0), Sackett's suggested values (middle panel; 1997), 
and Kerr & Lynden-Bell's recommendations (right panel; 1986). From bottom to top (A — B) increases from our best fitting value of 
(A - B) = 25.92 km s _1 kpc _1 , through Feast & Whitelock's (1997) best fit value (A - B) FW = 27.2 km s^kpc -1 ), to (A — B) = 
(A — B) pw + I-"" = 28.2 km s _1 kpc — 1 . For smaller values of (A — B ) and/or Ro, the rotation curve will decline even more with radius. 



The mass model fit to the rotation curve is shown as a solid line (§4.1). Also shown are the contributions from the individual components: 
the bulge (dotted line); the stellar disk (filled circles); the dark halo (long dashed line); and the gas layer (short dashed line). Notice 
that the contribution of the dark halo increases for smaller values of Ro and larger rotation sp eeds. The degree to which the disk is 
maximal (7; Oiling 1995) is indicated in each panel. Note that the total stellar matter (Sackett 1997) contributes somewhat more: ~74%, 
~81%, and ~81% for Ro = 7.1, 7.8, and 8.5 kpc, respectively. The non-linear behaviour of the rotation curve due to the gas in the Solar 
neighbourhood introduces a tiny wiggle in the total model rotation curves, similar effects are seen around R = O.6R0 and 1.4_Rq- 
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Ro values tend rise in the outer Galaxy. Figure [l] shows some 
typical models. 

A comparison of the Galactic rotation curve with ro- 
tation curves of external galaxies is best made when the 
radial extend is re-scal ed i n terms of the optical scale-length 
(hit = 2.5 ± 0.5 kpc; §4.1). The last measured point on the 
Milky Way's rotation curve lies at 2Ro, so that the extent 
falls somewhere between approximately 4.7/ir, (= 2 x 7.1/3) 
and 8.5/ia (= 2 x 8.5/2). It it not uncommon for galaxies 
to have rotation curves which decline somewhat between a 
peak at around 2.3/ir and 5/ir. This is in particularly com- 
mon in syste ms with large bulge-to-disk ratios^] like M31 
( |Kent 1989|), NGC 2683 & NG C 3521 dCascrtano fc van 



Gorko m 1991), and man y others (Rubin et al. 1985; Bege- 
man 1! '87 1_ Broeils 1992 ) . In terms of the universal rotation 
curve (Persic, Salucci & Stel 1996), 200 kms -1 is a divid- 



ing point: on average, systems which rotate faster decline 
beyond 2.3/ir, while for slower rotators the trend is to have 
a rising rotation curve. From the comparison with external 
galaxies we conclude that Oo = 180 ± 20 kms -1 , while a 
rotation speed at the solar circle as high as 220 km s _1 is 
highly unlikely (see also Merrifield 1992). 

The accuracy of the Milky Way's rotation curve - re- 
flected by the scatter of the data points in Fig. [l] - is much 
poorer than that routinely obtained for external galaxies. 
This is partly explained by the fact that for external galax- 
ies a large range in azimuth can be used, while for the Milky 
Way only two regions in Galactocentric longitude around 
I = 135° and I = 225° [cf. Eqn. @] are useful. Only in those 
regions is the iv-distance gradient large enough to rise above 
the noise induced by the tracer's random motions. Further- 
more, small errors in the tracers distances or small non- 
circular motions in regions away from I — 135° and I = 225° 
can result in large changes in the inferred rotational veloci- 
ties. Clear examples of this effect are the deviant points on 
our rotation curve between 1 and 1.2Ro which arise from 
r adial velocities of H2 regions between / = 90° and / = 120° 
(Brand & Blitz 1993). These da ta points are likely to b e 
affected by non-circular motions ( Dehnen fc Binney 1997]) 



induced by the nearby Perseus arm (Brand & Blitz 1993) 



and distance errors resulting from unaccounted for extinc- 
tion, or errors in the photometric parallaxes (Oudmaijer, 
private communications) . These large errors preclude the de- 
termination of the Oort constants from the Milky Way's ob- 
served rotation curve. As the best alternative, we will model 
the rotation curve and determine the Oort constants from 
the smooth model fits instead. Rather than using the cus- 
tomary linear or power-law fits to the rotation curve, we 
use the mass modeling procedure which has proven highly 
successful in external galaxies (§[l.l|). 



3 THE H I AND H 2 DISTRIBUTIONS 

In order to determine the radial variation of the column 
density of the molecular and atomic hydrogen, we need to 
be able to assign a distance to a parcel of gas observed at 



t The M ilky Way's bulge-to-disk ra tio is approximately 0.27, in 
K-band ( |piling fc Merrifield 1998a| ) 




Figure 2. The upper panel shows the face-on surface density 
distribution of the H I (dotted line), the H2 (dashed line), and 
the sum of two (full line) for the Milky Way. The lower panel 
shows the contribution these components make to the circular 
speed of the Galaxy, where a negative rotational velocity means 
that the radial force is directed outwards. For other values of Ro, 
the surface density is unaffected, but the rotation speeds vary 
approximately proportional to Ro. The radial force due to the 
ISM varies strongest in the regions around 0.6Rq, ^0; an d IARq. 



Galactic coordinates I and b, and radial velocity v r . Assum- 
ing that the gas is on a circular orbit, its distance d can be 
determined using the following relations: 



W{ 



Ro' 



R z 



sin I cos b 



e 



= Ro 2 + d 2 cos 2 b ■ 



2Ro d cos / cos b . 



(6) 



(7) 



While the radial velocity and W(-§-) can be determined ac- 



curat ely QMalhotra 1995| jMerrifie'ld 1992 



Brand & Blitz 



199S ) , the uncertainties in Ro and Oo lead to a significant 



uncertainty in the value of d. We have therefore performed 
this analysis using a range of values for Ro and ©0. Using 
the rotation curves which follow from the adopted values for 
Ro and Oo, we determined the volume and column densities 
of the molecular and atomic gas. For the inner Galaxy, the 
H I column d ensity was det ermined from the midplane vol- 
ume dens ity (Burton 1988 ) and the observed thickness of 
the layer ( Malhotra 199E| ).~The column densities for R > Ro 
were taken from Wouterloot et al. (1990), with 9.25 Mq pc -2 
at the Solar circle. The H2 column densities for the inner and 
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outer Galaxy were taken from Bronfman et al. (1988) and 
Wouterloot et al. (1990), respectively. The column density 
of H2 at the solar circle is 1.8 Mq pc -2 , and so, neglecting 
the column density due to the other phases of the ISM, we 
arrive at a total column density for the ISM at the solar cir- 
cle of 14.5 Mq pc~ 2 , whi ch includes 23.8% Helium by mass 
(Olive & Steigman 1995D - Note that while the column den- 
sity at fractional radius R/Rp is independent of the choice 
of Ro (Bronfman et al. 1988), the Milky Way's total gas 



mass does depend on this value. 

The H I and H2 column densities, and these mass com- 
ponents' contributions to the total rotation curve, are pre- 
sented in Fig. [| It is apparent that the sun resides in a part 
of the Milky Way where the gradient of the gas column den- 
sity is substantial, and that their corresponding contribu- 
tions to the gradient in the rotation curve are large. Clearly, 
this effect must be taken into account when interpreting the 
Oort functions of equations (|l|) and (^|). 

Note however that this effect is not limited to the Solar 
neighbourhood: the regions around 0.6i?o and IARq exhibit 
similar behavior. Examples of such non-monotonic rotation 
curves and H I surface dens ity distributions can be found 
in many external galaxies (Begeman 1987). Some galax- 
ies which have been observed with high linear res olution 



show strong radial H I density gradients: NGC 4244 ( Oiling 
1996a |)]NGC 891 (pccquacrt. Combes fc Viallefond 1997|) 



M31 (Prinks fc Burton 1984|), M33 ( peul fc van der Hulst 
1987|) , 4nd NGC 1560 QBroeils 1999b] ) all show similar radial 
density gradients to those found in the Milky Way. 



4 ROTATION CURVE, MASS MODELS AND 
THE OORT CONSTRAINTS 

In order to produce a complete mass model for the Milky 
Way, we must fit all the components - gas, stars and dark 
matter - to the observed rotation curve. 



4.1 Mass models 

In a previous section, we obtained the gas mass distribu- 
tion, so we now turn to the stars and dark matter. We 
have modeled the stellar component as the sum of a bulge 
and disk, whose light distribution is consistent with Kent's 
(1992) photometry, with the disk scalelengt h forced to lie in 
the range /ir=2.5 ± 0.5 kpc flSackett 1997|). The dark halo 



was modeled as a non-singular isothermal sphere. We refe r 



the reader to a related paper (Oiling & Merrifield 1998a) 
for a more detailed description of the mass models. 

We then solved for the free parameters in this model 
by fitting the observed rotation curve to the predictions of 
the mass model. This procedure that has been followed in 



many previous analyses (|Caldwcll & Ostriker 1981 




3ahcall, 


Schmidt & Soneira 1983; Ftohlfs & Kreitschmann 




Kuijkcn 


& Gilmore 1989|; Kent 1992[ 


Merrifield 1992; |Evans & Ji- 


jina 1994^ Alcock et al. 1995|; 


Gates, Gyuk & Turner 1995 


Dehnen & Binney 1997), but, to our knowledge, none of 



the previous studies have included the actual radial surface 
density distribution of the Milky Way's ISM. In addition, we 
imposed the constraint that the mass model had to repro- 
duce the observed stellar surface density of the Milky Way 



in the solar neighbourhood, E* = 35 ± 5 Mq pc 2 (Kuijken 
& Gilmore 1989). 

For each of the 50 i?o-©o combinations considered, we 
constructed nine mass models with /ir=2, 2.5, and 3 kpc, 
and three stellar mass-to-li ght ratios (T») such that E* is 
recovered to within ± 1-a (Oiling & Merrifield 1998a). We 
show nine of the ~450 mass models - chosen to be repre- 
sentative of the ensemble - in Figure |IJ 

4.2 Oort constants constraints 

From each mass model, we can calculate directly the Oort 
functions, A(R) and B(R), as defined by equations (jlj) and 
(^). Figure ^| shows these functions for models obtained from 
three combinations of Ro and Oo. It is clear that the Oort 
functions (and the various combinations thereof) vary sig- 
nificantly over even the limited range of radii shown. In the 
Solar neighbourhood about half of the variation in A(R) and 
B(R) can be attributed to the general Q(R)/R trend (dot- 
ted lines), while the local deviations are caused by the con- 
tribution of the gaseous component to the rotation curve. 
It is also apparent that the values of these functions vary 
significantly over the acceptable range of values for Ro and 

e . 

The functions also change to some extent when we 
adopt different acceptable combinations of the other param- 
eters in the mass model {ha, T», and the parameters of the 
dark halo). However, these modifications are all small. The 
error bar toward the right of each panel in Figure ^| indi- 
cates the full range of variations induced by such changes. 
Since the acceptable combinations of stellar and dark mat- 
ter components all reproduce the same large-scale smooth 
variation with radius of the rotation curve, the resulting 
Oort functions scarcely depend on the particular combina- 
tion adopted. The dependence on Ro and Qo arises because 
changing these quantities alters the largescale shape of Q(R) 
(see Fig. |l|). Thus, although the observed values for various 
Oort constraints will not distinguish between the details of 
the various degenerate mass models, they will allow us to 
decide which combinations of Ro and <do are consistent with 
the observed stellar kinematics. 

The observed values for the various Oort constraints, 
taken from Kerr & Lynden-Bell (1986; hereafter referred to 
as KLB86) and Hanson (1987; henceforth H87) are summa- 
rized in Table [j], and are also plotted in Figure ^| Since the 
Oort functions vary substantially with radius, it is vital that 
we compare the observed values to the functions over the ap- 
propriate range of radii. If we had not included the contribu- 
tion from the gas component (dotted lines in Figure ^), the 
Oort functions A(R) and B(R) would have changed steadily 
with Galactocentric radius [dA/dR ~ — 1.9,dB/dR ~ 2.5 
(km s _1 kpc _1 )/kpc]. However, the contribution of the gas 
flattens the A and B gradients substantially over the scale 
of the width of the H I crest which peaks at l.LRo. In this 
radial regime we can thus use the small distance approxima- 
tions - which require constant A and B values - for v r and 
fi, [Eqns. (§) and @]. 

Unfortunately, KLB86 obtained their estimates by av- 
eraging together a number of previous determinations, and 
it is not clear what range in radii the data for these vari- 
ous analyses came from. As KLB86 noted, it is even unclear 
if the individual estimates were all obtained from data at 
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Table 1. Th e Port Constants as determined b y several groups: KLB86 (Kerr & Lynden-Bcll 1986), PMB94 ( Pont, Mayor fc Burki 
1994^7^B93 (Brand & Blitz 1993), and FW97 (Feast & Whitclock 1997). We adopt the values as determined by H87 ( ^Tanson 1987 ). 
The units of A, B, (A — B), and (A + B) are km s -1 kpc - x , Rq is in kpc, and So in km s -1 . The ratio —B/(A — B) may be equated 
with the square of the tangential to radial stellar velocity dispersion (X 2 = o^/ctjj). Most stars of spectral type G through K fall in the 



"low X 2 " . M type stars fall mostly in the "high X 2 " category (see §4.2). 





A 


—B 


A — B = e /R 


A + B = (d@/dR)\ Ra 


-10B 
A-B 


Ro 


e 


KLB86 


14.5 ±1.3 


12.0 ±2.8 


26.4 ±1.6 


±2.5 ±3.1 


4.2 ±0.6 


8.5± 1.1 


220 ± 20 


BB93 


12.6 


13.2 


25.9 


-0.6 


5.1 


8.5 


220 


PMB94 


15.9 ±0.3 










8.1± 0.3 




FW97 


14.8 ±0.8 


12.4 ±0.6 


27.2 ± 1.0 


+2.4 ± 1.0 


4.5 ± 0.2 


8.5± 0.5 


231 ± 15 


H87 


11.3 ±1.1 


13.9 ±0.9 


25.2 ± 1.9 


-2.6 ± 1.4 


5.5 ± 0.5 






This work 


11.3 ±1.1 


13.9 ±0.9 


25.2 ± 1.9 


-2.6 ± 1.4 




7.1± 0.4 


184 ± 8 


low X 2 










3.5 ± 0.5 






high X 2 










5.0 ±0.5 







6 7 8 [kpc] 6 7 8 9 [kpc] 7 8 9 10 [kpc] 




Galactocentric radius Galactocentric radius Galactocentric radius 



Figure 3. The Oort functions, A(R) and B(R), as derived for three model rotation curves, together with the values of (A — B) and 
—B/(A — B). The units for A, B, and (A — B) are km s — 1 kpc -1 . The solid lines show the Oort functions as derived from the full 
mass model, while the dashed line shows what happens if the gas component is ignored. The error bar at the right of each panel shows 
the range of uncertainty due to variations in the other model parameters (T* and /ir). The various observational estimates of these 
quantities are also shown, with the horizontal error bars indicating the radial range over which the observations effectively averaged (see 
text for details). The upper three panels show values as determined by KLB86 (open squares) and H87 (open circles). The bottom panels 
show the values of —B/(A — B) as predicted by the velocity dispersions of "high X 2 " stars (M dwarfs and G-M giants; triangles) and 
"low X 2 " stars (G-K dwarfs; hexagons). 



comme: lsurate radii. We are therefore unable to estimate 



the appropriate radial range for these estimates. H87 used a 
more consistent set of data from the Lick Northern Proper 
Motion survey, consisting of 60,000 stars within ab out 1 kpc 
of the Sun. From their B — V colors [0.84 ± 0.5, ([Klemola 



Jones fc Hanson 1987 )] we infer that these stars are F2-K0 
dwarf stars ( jVlihalas fc Binney 198l| ) with a radial disper- 
sion of roughly 29 km s~ , so that the epicyclic radius for 
these stars is ~780 pc [cf. Eqn. (jB])]. The Galactic longitudes 
of these stars (I « 50 - 195 degrees; \b\ > 10°) indicate that 
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they are mostly located in the outer Galaxy: the horizontal 
bar on the H87 data points in Figure |§] represents the range 
of J? of a typical star in this sample [(Ro — 250, Ro + 1000) 
parsec]. Since this region roughly equals the area which is 
covered by the epicyclic motions of these F2-K0 dwarfs, they 
can be used to determine intrinsic Galactic properties on 
this scale. Furthermore, as discussed above, Hanson's use of 
the approximate fii relation [Eqn. (Q)] is justified due to the 
flatness of the A(R) and B(R) functions in this region. 

From Figure ra it is clear that the values of A and B 
derived by H87 are consistent with the Oort functions for 
small values of Rq and Go. Models without the contribution 
from the gas component yield larger values for A and B. ft 
is also noteworthy that the value of 2ARo based on deter- 
minations of W(R) in the inner Galaxy is 1.4 ± 0.25 times 



larger t han the value determined in the outer Galaxy (Mer- 



rifield -992); this finding is also entirely consistent with the 



Oort ^4-function curves shown in Figure |3j. 

An independent estimate for (A — B) comes from proper 
motion surveys in the directions / = 90 degrees and / = 
270 degrees. Objects in these directions lie at radii close to 
Ro out to quite large distances, and so no radial range is 
plotted for the KLB86 and H87 estimates for this quantity 
in Figure |^. In any case, the absence of a dB/dR term in 
this combination of Oort functions means that (^4 — B) is a 
smooth function of radius, so the importance of knowing the 
precise radial range of the determination is reduced. Once 
again, the H87 value for (A — B) is best matched by models 
with small values of Ro and Oo- 

Finally, we turn to the observed value of —B/(A — B). 
If the potential is azimuthally symmetric and the velocity 
dispersions are small, then this combination of the Oort pa- 
rameters is equal to the square of the ratio of the tangen- 
tial an d radial stellar velocity dispersions, (o^/or) 2 (Kui- 
jken & rremaine 1994 ). The a ccepted value for this ratio i s 



X 2 = (cr /cr R ) 2 =O.42 ± 0.06 ( [Kerr fc Lynden-Bell 1986[ ) 
However, upon close examination there appears to be a bi- 
modal distribution in X 2 : dwarfs of spectral type G-K have 
a low value [X 2 , — 0.36 ± 0.07), while M dwarfs and 
-M giant stars have a somewhat larRer dispersi on ratio 
{Xl M U 5 + 1 1 and X 2 r u = 5 + 07; f|riclhayf 



enceforth D65, and references therei 



3 



Mihalas & 



1965. 

Binney 1981)]. In fact, the large error on X^ M arises- as 
a resul t of the bimodality in the distribution of X 2 
Hawley 



Gizis 1995: 



X 2 ow — 0.36, jffejgh 



(Reid 



= 0.48), which 



is also clearly se en in the distribution of space motions of 
individ ual stars (Vyssotsky fc Janssen 1951 ; Dclhayc 1965 
and ref 
Binney 



j rences therein ;_ Delhaye 1965, Fig. 2 ; Mihalas k 
1981, Figs. 7.1 fc 7.3| ).' Furthermore, these figures 



show that the giant stars have a bimodal distribution as 
well. This result was confirmed by more recent investiga- 
tions (Mayor 1974; Oblak 1983), which found a minimum 
dispersion ratio for stars ~ 1 Gyrs old (X 2 ~ 0.36) while 
the dispersion ratios for older and younger stars are larger 
(X 2 ~ 0.46). Thus, rather than having to account for the 
average X 2 -value of 0.42 ± 0.06, we are faced with the task 
to explain the bimodal X 2 distribution. We estimate the lo- 
cation of the two peaks to lie at: X 2 — 0.35 ± 0.05, and 
X 2 = 0.5 ± 0.05. Figure H suggests that these results can 
be explained if the bimodality in the dispersion-ratio mea- 
surements is a result of the "dip" and the "hump" in the 
— B/(A — B) curve around 0.9i?o and 1.1-Ro, respectively. 
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Figure 4. The contours of equal likelihood as a function of 
Ro and 0o , calculated by comparing stellar kinematic constraints 
to the model-predicted Oort constants. The best-fit minimum-x 2 
values for Ro , ©Oi and their 1-tr errors are also indicated. 



In this scenario, the dwarf sample is dominated by stars 
originating within the Solar circle (small X 2 values) with a 
smaller number of dwarf showing the larger X 2 -value, while 
the giants used in the dispersion measurements come pref- 
erentially from the outer Galaxy. Since the stellar density 
increases towards smaller R with a scale-length that may 
be as small as 2 kpc, and since the observed distribu tion 
of giants is co ncentrated around (R — Ro) — 250 pc (Mc- 



Cuskey 1965), this may well be the case. The larger range 



in radii covered by the dwarf stars with their larger epicyclic 
excursions (indicated by the horizontal bars on the points 
in Fig. ^ also means that these stars are more influenced 
by the dip in the —B/(A — B) curve at R ~ 0.9i?o, ex- 
plaining the lower observed average value. Once again, only 
models with low values of Ro and Oo are consistent with the 
observations. 



5 DISCUSSION 

It is apparent from the above analysis that the local to- 
pography of the rotation curve induced by the small-scale 
structure in the gas component has a significant impact on 
the interpretation of the stellar-kinematic Oort constraints. 
Once this phenomenon is taken into account, it is possible 
to obtain model rotation curves for the Milky Way that are 
consistent with essentially all the Oort constraints. 

It is also clear from Figure ^ that a consistent picture 
only emerges for certain values of Ro and ©o. We can formal- 
ize the constraints that this analysis puts on the values of the 
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Galactic constants by calculating a x 2 statistic comparing 
the observed combinations of the Oort constants in Figure H 
to the values predicted by the models. Since this comparison 
requires that we know the radial range over which the ob- 
servations have been made, we use the observations plotted 
in Figure ^ for which we have an estimate of this range (see 
above), together with the H87 estimate for (A — B). 

The \ 2 statistic comparing these five observations to 
suitably radially-averaged models were then calculated for 
a range of values for Ro and Oo- The best-fit (minimum- 
X 2 ) values for these parameters are Ro = 7.1 ± 0.4 kpc, and 
Oo = 184 ± 8 km s _1 . The complete two-dimensional confi- 
dence regions are plotted in Figure The confidence levels 
measure the probability that any given values for Ro and So 
are consistent with the observed Oort constraints. Thus, for 
example, the official IAU-sanctioned values of Ro = 8.5 kpc 
and Oo = 220 kms _1 are ruled out at the 99% confidence 
level. It is, however, notable that a value for Ro as low as 7.1 
kpc is entirely consistent with the only direct distance deter- 
mination to the Galactic center em ploying prop er motions 
of H 2 masers, Ro = 7.2 ± 0.7 kpc flReid 199: 

The proper mot i on of Sgr A 
1987; E acker 1996; Eteid 199£), has been used to used to 
constrain (A — B) as well ( |Gould fc Ramirez 1997 ). These 
quantities are related via the formula, 



HSgrA* (Backer & Sramek 



HSgrA* = ' 



A-B 



-J milharcsec 



year 



(8) 



4.74 4.74.Ro 

where vo is the Solar motion relative to the local standard 
of rest (vo = 5.3 ± 1.7 kms -1 ; Binney et al. 1997) and 
vsgrA* is the velocity of Sgr A* with respect to the Galactic 
centre (vs gr A' < 11 km s _1 ; Gould & Ramirez 1997). Using 
the range of acceptable values for (A ~ B) derived from 
our mass models, we find that the proper motion of Sgr A* 
should lie in the range fisgrA* = —5.95 ± 0.36 milliarcsec 
per year, which is entirely consistent with the observational 
constraints. 



5.1 Other determinations of the Galactic 
constants 

Inspection of Table |l| reveals that the Oort constants as 
derived by different authors can vary by up to 4.2-er. We 
have used the determination of the Oort constants which are 
based on ground based proper motions of relatively nearby 
stars (d < 1 kpc; H87) in the outer Galaxy. Cepheids probe 
a much larger range in distance and can be used to deter- 
mine the O ort constants as well, but yield significantly differ- 
ent values flPont, Mayor fc Burki 1994, henceforth PMB94j ; 



Feast fc Whitelock 1997, hereafter refered to as FW97) 



Let us first examine the region within 1 kiloparsec from 
the Sun where the non-linearity of the Oort functions is 
most prominent. As discussed above, the H87 determined A 
and B in the "flat" region extending from 0.957?o to 1.157?o: 
Ah87 = 11-3 ± 1.1 km s _1 kpc _1 . The Cepheids on the other 
hand are located preferentially in the inner Galaxy (PMB94, 
their Fig. 4) where our models yield an average A of 15.4 ± 2, 
entirely consistent with their value of 15.9 ± 0.3. From Fig- 
ure [] it is also clear that beyond 1 kpc the non-linearities in 
A(R) and B(R) are much smaller. In order to assess whether 
the PMB94 & FW97 results are consistent with our models 
at large distances, we have compared the radial velocities 



and proper motions as predicted by the PMB94 & FW97 
models with our models. Since the PMB94 & FW97 models 
are fitted directly to their data, we are in effect comparing 
our models with their data. Our model rotation curves yield 
exact radial velocities with equation (^), and exact proper 
motions using the formula 



/ Q(R) (Ro cos / - d) cos b n , N 
Hi = I — ^ Oo cos I 



4.74d 



(9) 



in units of mi lli arcsec per year, with the distances in units 
of kiloparsec (Mihalas & Binney 1981). In the upper panels 
of Fig. ^ we present the proper motion curves as modeled 
by FW97 (symbols) and some of our models (lines) for a 
range of values for Ro, Oo, and d. In the bottom panels 
we plot the difference between the two. Our model which 
best fits the Oort constant constraints (i?o=7.1 kpc, ©o=184 
km s _1 ) is displayed in the leftmost column; the agreement 
between the FW97 proper motions and our predictions is 
excellent. The worst agreement between data and model is 
found for large distances towards the Galactic centre, where 
the rotation curve is rather ill-defined as a result of bar- 
induced non-circular motions. Furthermore, we note that 
models with larger and smaller values for Ro fit the FW97 
data progressively worse. We draw the same conclusion from 
a similar comparison of the PMB94 radial velocity data with 
our models. 

From Figure [s| we can also see that the vertical scatter 
seen in observational I — /i; plots is not due to some intrinsic 
scatter, but rather results from the fact that the data oc- 
cupy a finite range in radii. Clearly the approach followed 
by PMB94 & FW97 - fitting model v r and hi curves which 
include distance effects - is superior to the application of 
the approximate equations [Eqns. (^) and (^)]. 

Given that the best-fit values of the Galactic & Oort 
constants in this analysis are so different from the Cepheid 
derived values it behooves us to explain why our best model 
nevertheless is able to fit the radial velocity and proper mo- 
tion data (Fig. ^) so well. The reason for this is that the 
Milky Way's rotation curve and the Galactic & Oort con- 
stants are secondary parameters: a wide range in secondary 
parameters is consistent with the primary observable W(-g^) 

(Fig. [j]). So assuming that the W(-j^) curves as measured by 
Brand & Blitz (1993), Merrifield (1992), and Pont, Mayor & 
Burki (1994) are a fair sample of the intrinsic Galactic ve- 
locity field, it is not surprising that our smooth model fits to 
the first two data sets are consistent with PMB94's Cepheids 
data. Varying the Galactic constants has only a small effect 
on the proper motions, as can be seen by rewriting equation 
(^) in terms of the primary observable W 



Hi - 




Ro 



(A-B) 



(10) 



Changing (A — B) shifts the /i; curves up and down, but the 
observed error on (A — B) corresponds to a small vertical 
shift of Hi, which is only significant in the Galactic centre 
and anti-centre regions. Likewise, changing Ro affects the 
ratio - and hence W - only slightly [through Eqn. (Q)]. 
Like a change in (A—B), adjusting Ro shifts the Hi curves up 
and down. The only longitudinal dependence comes from the 
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Figure 5. Proper motions (top row) as calculated using the Oort and Galactic constants derived by FW97 (Ro = 8.5 kpc, A = 
14.82,5 = —12.37 km s — 1 , and d 2 Q/dR 2 = d 3 @/dR 3 = 0), and some of our models. In the lower panels we display the difference 
between the proper motions predicted by our models and those predicted by FW97. The units for the vertical axes are km s _1 kpc — 1 . 
All curves were calculated for 6 = 0. The error bar in the bottom panel equals the quadratic sum of the errors of FW97's determinations 
of A and B. In the leftmost panel we present our model which best fits the Oort constants constraints (Ro =7.1 kpc, ©o = 184 km 
s _1 ), while the other two panels were calculated for the Sackett (1997; middle panels) and Kerr &: Lynden-Bell (1986; right columns) 
recommendations for the Galactic constants. Again, worse correspondence between models and "data" occurs if (|i?o — T. 1 1 ) increases. 
In our models, assuming a different value for Oq will shift the m curves slightly up or down by A@q/R. 



first term, which can be large for small distances d. Thus, the 
differences between the proper motion curves for different 
choices of Ro & Go in Fig. [] are not very large. However, the 
residuals do increase if either Ro or Go or both are adjusted: 
our preferred model fits the proper motion data best. 



5.2 Uncertainties in the W(R) curve 

The derived form for the Galactic rotation curve depends 
on the adopted W(R) curve. In this paper, we have based 
our analysis on the radial velocities of H II regions (Brand 
& Blitz 1993) and the H I W(R) curve (Merrifield 1992). 
However, one could also infer this curve from radial veloc- 
ities of Cepheids. As Dehnen & Binney (1997) have noted, 
the W(R) curve derived from Cepheids drops significantly 
more rapidly with R than that obtained from the data sets 
exploited here. For such a rapidly-falling form for W(R) to 
produce a reasonably flat rotation curve, one must adopt a 
large value of Go - It is, therefore, unsurprising that Metzger, 
Caldwell & Schechter's (1998) analysis of Cepheid kinemat- 
ics results in models of the Milky Way with Go as high as 
240 ± 10km s" 1 . Their modeling obtained a value for Ro of 
7.7 ± 0.3 kpc; the corresponding value for (A — B) — G0/-R0 
is 31±2km s _1 kpc -1 , which is inconsistent with the various 
stellar kinematic constraints summarized in Table |l|. 

Clearly, these internal inconsistencies mean that there 



must be systematic errors in at least some of the available 
data. The unknown nature of these errors means that there 
can be no objective criterion for choosing between the var- 
ious methods. The only way that a final answer can be 
reached will be from the analyses of other kinematic trac- 
ers. Through such studies, one might hope to identify and 
prune the discrepant data sets, obtaining a consistent pic- 
ture from the remaining majority. With only three "reliable" 
techniques for obtaining W(R), we are not yet in a position 
to carry out such a democratic procedure. 



5.3 The effects of non-circular motions 

One potential source of error in this analysis comes from 
our assumption that the orbital structure of the Milky Way 
has an azimuthal symmetry. The presence of spiral struc- 
ture means that this assumption is almost certainly an over- 
simplification. Non-circular motions will have two effects on 
the above analysis. 

First, in Section^ we used kinematic distances to derive 
the distribution of gas. These distance assume that the gas 
follows circular orbits, and so the inferred distribution will 
be distorted if there is significant non-circular streaming. 
However, dynamical models of spiral structure show that 
the locations of high density ridges are quite accurately re- 
produced by the naive kinematic modelling, and so the ra- 
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dial positions of the density peaks and troughs apparent in 
Figure ^ should be approximately correct, and the conclu- 
sion that we are located in a region where this density varies 
rapidly with radius still holds. Second, the non-circular mo- 
tions hamper the determination of the intrinsic W(R) curve 
from which the rotation curve, the Oort constants, and the 
Milky Way's mass models are derived. With sufficiently large 
azimuthal coverage and dense sampling, such regions can 
be identifi ed [e.g., the Perseus arm between I — 90° and 
I = 120°, ( {grand fc Blitz 1993b ], and discarded. However 
note that when streaming motions have maximal effect on 
the observed radial velocities, the effects on the proper mo- 
tions will be minimal, and vice-versa. Thus, by combining 
the information from the line of sight and transverse mo- 
tions, all azimuths can be properly sampled, which should 



lead to a significantly more accurate determination of the 
Milky Way's rotation curve. One could then use the devia- 
tions from the general axisymmetric velocity field to study 
the non-circular motions in great detail and infer the mass 
associated with the spiral arms. 

We have seen that the distribution of the ISM has a 
large effect on the determination of the Oort constants in 
the Solar neighbourhood. However, the non-circular stream- 
ing motions have not been mapped out within one kpc from 
the Sun because random motions dominate the radial ve- 
locity and proper motion measurements of individual tracer 
objects. Although Hanson (1987) overcame this problem by 
using over 60,000 stars in his determination of the Oort 
constants, his data set cannot be used to investigate non- 
circular motions in the Solar neighbourhood as no distances 
to individual objects are available for this survey. 

We should keep in mind that if streaming motions are 
present to a significant degree, the stellar orbits will become 
more complex, and so the Oort constraints derived from 
simple stellar kinematic measurements will also be com- 
promised. Clearly, in order to allow for these phenomena 
fully, one must compare a complete dynamical model for 
the Milky Way to the observational data. However, appli- 
cation of Occam's razor to this problem suggests that the 
first step should be the comparison of observations to the 
simpler azimuthally-symmetric model. The fact that we can 
satisfactorily describe the Oort constraints with this model 
implies that there is no direct evidence for non-circular mo- 
tions in Hanson's Solar neighbourhood data. If, however, 
future observations establish that the true values for the 
Galactic constants lie outside the acceptable region in Fig- 
ure ^, then at least one of the assumptions underlying our 
models must be wrong. For example, it might be that signif- 
icant non-circular stellar motions are present. Alternatively, 
the mass distributions of the stars, or indeed the dark mat- 
ter may be non-smooth and contribute significantly to local 
gradients in the rotation curve. 



6 SUMMARY 

We have used distance and radial velocity data available in 
the literature to investigate the rotation curve of the Milky 
Way. The accepted uncertainties in the Galactic constants 
of ~10% allow for a wide range in rotation curve shapes. We 
have fitted mass models to about 450 rotation curves which 
are consistent with the radial velocity data and determined 



the Oort functions A(R) and B(R) for these model rota- 
tion curves. Typically, the Oort functions vary at a rate of 
a few km s -1 kpc -2 . In several locations in the Milky Way, 
the more objective model rotation curves show sufficient fine 
structure to affect the Oort functions significantly. The ring 
of H I just beyond the Solar circle is massive enough to flat- 
ten A(R) and B(R) in the first kiloparsec beyond the Solar 
circle. Comparison between the model A(R) and B(R) func- 
tions and Hanson's (1987) determinations of the Oort con- 
stants places tight limits on the distance to the Galactic cen- 
ter and the Milky Way's rotation speed at the position of the 
Sun: Ro = 7.1 ± 0.4 kpc, and 6 = 184 ± 8 km s" 1 . These 
results are consistent with the only direct determination of 
the distance to the Galactic center, _Ro = 7.2 ± .7 kpc, by 
Reid (1993), and the proper motion o f Sgr A* (Backer & 



Sramek 1987; Backer 1996: Reid 199£). With these Galac- 



tic constants we find that the rotation curve of the Milky 
Way declines slowly in the outer Galaxy. At larger distances 
from the Sun these models are entirely consistent with re- 
cent determinations of the Oort constants based on radial 
velocities and proper motions of Cepheids. 
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